Differential decay rate for B irliy semileptonic decays* 

JLQCD Collaboration: S. Aoki^ R. Burkhalter'^-'', M. Fukugita^ S. Hashimoto'^, K-I. Ishikawa'^, 
N. Ishizuka''^^, Y. Iwasaki^*-'^, K. Kanaya'^''", T. Kaneko'^, Y. Kuramashi'^ , M. Okawa'^, T. Onogi°, 
S. Tominaga'^, N. Tsutsui'^, A. Ukawa''^^, N. Yamada'^, T. Yoshie ^-'^ 

^ Institute of Physics, University of Tsukuba, Ibaraki 305-8571, Japan 

^ Center for Computational Physics, University of Tsukuba, Tsukuba, Ibaraki 305-8577, Japan 
Institute for Cosmic Ray Research, University of Tokyo, Kashiwa 277-8582, Japan 
High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 

° Department of Physics, Hiroshima University, Higashi-Hiroshima, Hiroshima 739-8526, Japan 

We present our study on B nli/ semileptonic decay form factors using NRQCD action for heavy quark. In 
the analysis, we use the form factors /i (v-k-^) and /2 (v-k-^) defined in the context of heavy quark effective theory 
by Burdman et al.. Since the NRQCD action for heavy quark respects the heavy quark symmetry, our results 
are described by the HQET form factors most naturally. From a quenched lattice QCD simulation at /3=5.9 on a 
16^ X 48 lattice, we obtain the form factors fi{v ■ fc^) and f2{v ■ /c,r), and find that their l/rriB correction is small. 
The limit of physical heavy and light quark masses can be reached without introducing any model function, and 
we obtain a prediction for the differential decay rate dF/dq^. We also discuss the soft pion limit of the form 
factors. 



1. Introduction 

Model-independent calculation of the B meson 
semileptonic decay form factors is essential for ex- 
tracting \Vub\, which is one of the most poorly 
determined elements of the CKM matrix. While 
the experimental results are being improved by 
the high statistics data from the B factories, there 
are also theoretical progresses in lattice QCDI- 
^ , which has a potential of making a precise pre- 
diction of the form factors near zero recoil. 

The success of lattice QCD for heavy-light 
physics especially in the decay constant may be 
attributed to the use of the heavy quark (or non- 
relativistic) effective theory (HQET). The effec- 
tive lattice actions constructed as a systematic 
expansion in 1 /mg provide a control over the sys- 
tematic errors from the large heavy quark mass. 
Furthermore, by identifying a quantity which has 
a well-defined heavy quark mass limit using the 
HQET, one may consider an expansion in 1/tob 
and study the effect of finite heavy quark mass 
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in a systematic way. Such a quantity for the de- 
cay constant is /s^/ms, for which many lattice 
calculations found that the 1/tob correction is 
significant. 

In this study we take the same strategy for B — > 
nlv semileptonic decay and use HQET motivated 
form factors /i(u ■ k^^) and f2{v ■ k^^) proposed by 
Burdman et al. as ^] 

-^{n{k.)\V^^\BipB)) 

k^^ 

- h{v■k^)v^^ + h{v■k^)^, (1) 

where and pb are four-momenta of tt and B 
mesons in the external states, and = Pb/''^b 
is a four-velocity of the B meson. The form fac- 
tors /i(w ■ kjr) and f2{v ■ k-^) are well-defined in 
the niB —foo limit, and it is natural to consider a 
l/niB expansion around that limit. We also note 
that fiiv-k^y) and /2(w'fc7r) are implicit functions 
of light quark mass rrig. The relation to the con- 
ventional form factors f^{q^) and f^{q^) is given 
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f2{v ■ k^) ^ fi(v ■ K) 



V ■ fcjr 



niB 




niB 



[h{v-K) + klh{v-K) 



(2) 



fi{v ■ k^) + /2(w • k^) 



(3) 



where we define kf^ ~ fc^/(u • fc„) and = m\ 
ml - 2mB{v ■ fc^). 
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2. Lattice calculations 

We have performed quenched simulations at 
(3=5.9 on a 16^ x 48 lattice, for which the lattice 
scale determined from the string tension is 1/a = 
1.64 GeV. We use the NRQCD action including 
0{\/mQ) terms for heavy quark and the 0{a)- 
improved light quark action with csw calculated 
at one- loop csw = 1-58. For the heavy quark 
mass we take oMq = 5.0, 3.0, 2.1 and 1.3, which 
cover the h quark mass. The light quark mass 
parameters, and some of the heavy quark mass, 
are the same as in our previous study for Jb ||]- 
We accumulated 2150 independent gauge config- 
urations to reduce statistical error for the signals 
with finite spatial momenta. Even with this large 
number of statistics, the signal for heaviest heavy 
or lightest light quark is not clean enough to ex- 
tract ground state. 

The matching of the heavy-light vector cur- 
rent = q^i^Q is done at one-loop by Morn- 
ingstar and Shigeniitsu |^ and by Ishikawa 
et al. M, in which the lattice opera- 
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components. In this work, we use the V-scheme 

coupling av{q*) with q* = 1/a for the coupling 

constant. 



Figure 1. The HQET motivated form factors in 
unit of a^/^ for four different heavy quark mass. 



In the measurement of the matrix element we 
take several combinations of initial and final me- 
son momenta. The maximum momentum we can 
measure is (1,0,0) for both initial and final mo- 
menta in the unit of 27r/(16a). As a result we can 
roughly cover a region 0.5~0.7 for av ■ k^^, where 
the lower limit is given by a pion mass am^. By 
taking an inner product with of both sides of 
(Q) , we extract a linear combination of form fac- 
tors fi{v ■ /cjr) + /2(f • krr). Another convenient 
basis is f2{v ■ k^r), which is obtained by taking a 
inner product with an unit vector perpendicular 
to v'^. 



3. HQET form factor results 

For the HQET motivated definition of the form 
factors it is natural to expand the heavy 
quark mass dependence in terms of I/ttib- The 
light quark mass dependence can be expressed by 
a Taylor expansion in mq when v ■ k^^ is fixed. 
Thus, we take the following form to fit the ma- 
trix elements 



fi{v ■ k„) + f2{v ■ k^) 

= Ao + mqAi + —A2, 
niB 



(4) 
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Figure 2. The HQET motivated form factors in 
unit of a^/^ for four different light quark mass. 



f2{v-K) = Bn + mgBi + B2, (5) 

where the expansion coefficients Ai and Bi are 
functions of u • fc,r- The expansions are trun- 
cated at first order, since we find only mild mass 
dependences, as we shall see. In practice, we 
parametrize v ■ k^r dependence of the parameters 
by a Taylor expansion around v-k^,- — Eq to a. cer- 
tain order. We take qEq to be 0.5, and expand 
Aq through quadratic term, and Ai, Bq through 
linear term. Others are taken to be constant. 

In Figure |l] the form factors /i (w • fc^ ) -I- /2 (w • fc^ ) 
and f2{v ■ fc^) are plotted for four different val- 
ues of aMo- We find that fi{v ■ k.^) + f2{v ■ fc^r) 
(filled symbols) has a significant negative slope 
in • fc^, while f2{v ■ fc^) (open symbols) is con- 
sistent with constant. Dotted curves in the plot 
represent the fit (^ and (H) . It is interesting that 
there is almost no heavy quark mass dependence, 
and therefore the 1 /to_b correction to these form 
factors is very small in contrast to the decay con- 
stant, for which a large slope in l/ms was found 
in lattice calculations. 

At a fixed heavy quark mass (aMo = 1.3), the 
form factors are plotted for different values of 
light quark mass in Figure |[ where the fit curves 
are also shown by dotted curves. We observe a 



Figure 3. Differential decay rate for B ttIv. 



mild dependence on the light quark mass. For 
/i(u • fc^) -I- /2(u • fcjr), a large fraction of the de- 
pendence comes from a trivial shift in u-fc7r, which 
is E'tt ~ -y/mir+kj in the B meson rest frame. 
An extrapolation to the physical or chiral light 
quark limit using the fit formula (^ and (|^) is 
shown by dashed and solid curves respectively. 

4. Differential decay rate 

From the fit results for the HQET form factors 
we may obtain the usual form factors f^{q^) and 
f^{q^) using (||) and (^. The region where our 
lattice data is available is 18-22 GeV^. The differ- 
ential decay rate dV /dq^ is then obtained through 



dq^ 



247r3 



|K.n/+(g^)P 



(6) 



which is plotted in Figure ^. We note that only 
the results in the region q^ = 18-22 GeV^, which 
is given by filled symbols, are obtained by inter- 
polating the simulation data in v ■ k^^, while the 
results outside the region involve an extrapola- 
tion of Ai and Bi in ?; • fc,r. Therefore, we should 
always keep in mind that the latter are subject to 
systematic errors from the choice of fitting func- 
tions. 

The error bars in the plot show the statistical 
error only. Systematic errors, such as the dis- 
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Figure 4. Soft pion limit of foiq'^ax) = 

[fi{v-kTr) + f2{v-kTr)]- The dashed line is a 
linear fit in {arriTr)'^, while the solid curve includes 
the term (amTr). A result of the fit is given by 
an open triangle, which should be equal to Jb //tt 
(filled triangle) in the soft pion theorem. 



cretization effect 0((ak7r)'^) or the pcrturbative 
error 0{a^), are under investigation. 

5. Form factors in the chiral limit 

Despite the uncertainties in the extrapolation 
ti • fcjr — > 0, it is important to study the form 
factor in the soft pion limit (mjr, A:^ ^0), since 
the HQET and chiral symmetry predict a relation 



/l(0)+./2(0) = 



max / 



(7) 



In 



where fs is the heavy-light decay constant. 
Figure ^ we compare the result of the fit (^) 
shown by an open triangle with the lattice cal- 
culation of /b/Ztt (filled triangle) Since the 
systematic errors, such as the uncertainty from 
V ■ fc^r extrapolation and higher order pcrturba- 
tive corrections, are not yet included in the er- 
ror bar, we consider that they are in reasonable 
agreement. We also plot two extrapolations of 

foilmax) = 2/VTOi'[/l(u • fc^) -f /2(w • fc^)] with 

the usual linear form in (am^)^ (dashed line) and 
with a quadratic fit in (am^) (solid curve). Since 
the value of w • fc^ varies in the extrapolation, lin- 



ear dependence in appears implicitly and the 
usual linear extrapolation in (or mq) is no 
longer justified Although the effect of the lin- 
ear term in am^ is very small and is only seen 
at the lightest quark mass, it raises the soft pion 
limit for the quadratic fit and makes it consistent 
with the result from extrapolation using the fit 

i). 

Near the zero recoil limit, the HQET predicts 
the B* pole dominance |^ 

f2{v ■ K) 



V ■ fc^ 



^ gBB-TT — TT? ; — ; ^ l°J 

S/tt V ■ k^, + niB' - mB 

where gBB*-n denotes the B*Btt coupling in the 
heavy-light meson chiral effective theory. Since 
the hyperfine splitting tub' — mB is small com- 
pared to w • /ctt , we can approximate its functional 
form by a constant in our data region. Then, we 
obtain gBB'irifB* ^/aniB' /2f-^) = 0.32(18), which 
gives gBB'-K = 0.27(15). It agrees with the phe- 
nomenological value extracted from D* — > Dtt 
decay 0.27(6) §. 
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